Dilation-induced phases of gases absorbed within a bundle of carbon nanotubes 
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A study is presented of the effects of gas (especially H2) absorption within the interstitial channels 
of a bundle of carbon nanotubes. The ground state of the system is determined by minimizing the 
total energy, which includes the molecules' interaction with the tubes, the inter-tube interaction, and 
the molecules' mutual interaction (which is screened by the tubes). The consequences of swelling 
include a significant increase in the gas uptake and a 3 per cent increase in the tubes' breathing 
mode frequency. 



PACS numbers:68.45.Da, 61.48+c, 82.65.My 

Considerable attention has focused in recent years on 
the absorption of H2 (and other gases) in various forms 
of carbon, especially high surface area materials, because 
of the potential such materials present for efficient stor- 
age, isotope separation and other applications. Numer- 
ous techniques of interfacial science are being applied to 
study the remarkable properties of this system, fn the 
case of carbon nanotubes (NT), the unusual geometry 
presents the possibility of novel phase transitions. The- 
oretical work has explored a variety of transitions and 
proposed ways to observe them experimentally . 

In this paper we propose such an unusual transition for 
hydrogen and other gases, associated with the dilation of 
a lattice of nanotubes. The most striking feature of a sin- 
gle nanotube is its quasi-one dimensional (ID) character, 
which makes contact with a growing body of ID theory 
Q. This geometry stimulated Gordillo, Boronat and Ca- 
suUeras Q to study a system of H2 molecules in ID and 
quasi-lD (in which case the molecules' small amplitude 
motion perpendicular to the tube's axial (z) direction was 
taken into account). They found density-dependent tran- 
sitions at temperature T = to a liquid phase and then 
to a novel ID solid phase. These results are qualitatively 
similar to those found by Boninsegni and coworkers [QJ^] 
for ^He in both ID and within an ordered lattice (bun- 
dle) of parallel NTs. If we consider such a bundle, then 
the interactions between molecules in adjacent intersti- 
tial channels (IC) allows a fully 3D transition to occur. 
However, the large spacing (^ 10 A) between adjacent 
channels means that the inter-IC interaction is weak and 
the predicted ordering temperature Tc is somewhat lower 
than the energy scale set by the well depth e of the pair 
potential The transition is one to a very anisotropic 
liquid (in the case when the molecules are assumed to 
move freely along the channels) or crystal, depending on 
the system 

Two new aspects of the absorption problem alter this 
situation both qualitatively and quantitatively. One is 
that a study of three-body interactions found 0] that an 
effect of the NTs' dynamic polarization is to significantly 
screen the interaction between molecules, with drastic 



consequences If one includes only intra-IC interac- 
tions, the ID II2 ground state in this environment is a 
very low density gas instead of the high density, strongly 
self-bound liquid mentioned above. However, inclusion 
of the inter-IC interactions does yield a weakly bound 
condensed state [||. The other new ingredient in this 
problem is some suggestion and tentative evidence , 
that adsorbed gases cause the bundle of NTs to swell. In 
this paper, we consider the problem of such a dilation 
of the lattice of NTs. We report four predictions associ- 
ated with this dilation: a greatly increased binding of the 
hydrogen, a significantly higher critical temperature for 
the condensed state, a measurably larger lattice constant 
of the NT array, and a higher breathing mode frequency 
for the tubes. Several other consequences merit future 
study. Especially interesting is the effect of this dilation 
on the electrical transport properties of the NTs, which 
have been found to be significantly altered by adsorption 
While we focus on the case of H2, we have found 
and summarize below qualitatively similar conclusions 
for the adsorption of He, Ne,Ar, and CH4. The unifying 
concept is that lattice dilation permits small molecules to 
increase their IC binding energy significantly without a 
substantial increase in the inter-tube interaction energy. 

The logic of our calculation is the following. Consider 
an infinite array of NT's, each of which has radius R 
and a large length L, and whose parallel axes intersect an 
orthogonal plane to form a triangular lattice of spacing 
rfo = 17 A at T = 0; we are ultimately interested in 
the thermodynamic limit of infinite L. We assume that 
between the tubes, within the IC's, there exists a ID 
density p of H2 molecules, so that on average there are 
N = pL molecules per channel. We then minimize the 
total energy of the system by allowing the lattice to dilate 
to a separation d. Such a resulting situation turns out to 
be stable only if p exceeds a threshold density, pc- This 
threshold density and attendant dilated lattice represent 
the ground state of the system. The finding that any 
density lower than pc is impossible is analogous to the 
familiar absence of a stable, uniform low density regime 
for conventional systems in free space at T = 0, i.e., the 
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ground state of all 3D systems is a condensed liquid or 
solid. As in the latter case, at finite T a low density 
gas phase does exist within the NT bundle, as does a 
regime of coexistence between this dilute gas and the 
dense liquid. 

The calculations themselves are straightforward, 
thanks to several simplifying assumptions which should 
not greatly affect the major conclusions; these include 
the neglect of any z dependence in the potential expe- 
rienced by the H2 molecules. The total energy of the 
system, per unit length, per IC, indicated by e, is writ- 
ten in terms of the single molecule's energy per particle 
(ei), of the H2-H2 interaction energy per particle (ei„t), 
and of the elastic energy per unit length representing the 
interaction between pairs of adjacent NTs as: 

e = E/L = p[ei{d) + e„t(p, d)] -f ifc(d - dof (1) 

The coefficient k — 1740 K A^'^ is derived from the 
semiempirical NT interaction constant of Mizel et al . 
This value is consistent with the experimental value of 
the interlayer force constant and interaction energy of 
graphite. The fimction eint{p, d) is our variational upper 
bound to the ground state energy per particle of fully 
interacting H2, computed with the screened interaction. 
This quantity is the sum of the "exact" ground state 
energy of ID H2 (computed by Boronat and Gordillo 
, using the diffusion Monte Carlo method) and the 
small intcrchanncl interaction, computed in the Hartree 
approximation. Finally, £i{d) is the ground state energy 
per H2 molecule subjected to the potential energy within 
an IC of a bundle of tubes which are spaced a distance d 
apart. We deduce this potential energy from the model of 
Stan and coworkers and solve the Schrodinger equa- 
tion using the diffusion method in order to obtain ei(c?). 
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FIG. 1. Azimuthally averaged potential energy and 
ground state eigenvalue (bottom panel) and probability den- 
sity (top panel) of H2 as a function of perpendicular distance 
from the center of the interstitial channel. Dashed curves 
correspond to the undilated NT lattice and full curves cor- 
respond to the lattice dilated by 1 per cent (the predicted 
ground state of the system). 
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FIG. 2. Chemical potential as a function of the ID den- 
sity p of H2 (full curve). The horizontal line exhibits the 
equilibrium (two-phase) behavior obtained from Maxwell con- 
struction which determines the uptake threshold density pc. 
Dash-dot curve is hypothetical result in the absence of dila- 
tion. 

Figure |l| shows the potential energy and wave function 
for the cases d = do = 17 A and a slightly larger value, 
d = 17.18 A. A key difference between these two cases is 
that the ground state of the dilated lattice has a much 
lower value of ei (d) . The reasons for this decrease are the 
lower potential energy and zero-point energy (ZPE) as- 
sociated with the increased distance of the H2 molecules 
from their neighboring NT's. This decrease in ei (equiva- 
lently, the quantum pressure, which in H2 is near 100 atm 
at threshold) is what drives the NT's apart |l6|. Figure 
§ shows the chemical potential (m = = f|) of the hy- 
drogen as a function of p, taking into account the lattice 
dilation. 



-100 

z 

-400 



-700 




0.1 0.2 PcO.3 



2 



FIG. 3. Energies as a function of the ID density p of H2. 
Filled circle is the NT-NT elastic interaction energy, dash-dot 
curve is the H2-H2 interaction energy per molecule and dashed 
curve represents the H2 molecule's eigenvalue in the dilated 
lattice. Full curve is the sum of the above contributions. No- 
tice that its minimum occurs at pc. 

Also shown is the analogous result in the case of an 
undilated lattice. In the latter case, the curve exhibits 
an extremely shallow minimum near p = 0.1 A~^, which 
is not discernible in the figure; the binding energy in 
that case is ~ O.IK In the present case of dilation, 
the chemical potential is seen to decrease rapidly with 
increasing p because the expanding lattice becomes pro- 
gressively more attractive, as seen in Fig. |l|; the eigen- 
value of this potential has a minimum at lattice constant 
17.55 A, corresponding to a very large density. However, 
the mutual II2 repulsion at such high density is huge. The 
ground state energy of the system is obtained by minimiz- 
ing the grand potential, uj — e — fip. This is equivalent to 
a Maxwell construction in the chemical potential-density 
plane, as seen in the figure. It corresponds cquivalently 
to setting the ID pressure of the H2 equal to the stress 
of the tubes, associated with the change of their total 
interaction energy with respect to density. 
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FIG. 4. Lattice constant of the NT array as a function of 
ID density for ''He, H2 and D2. The discontinuities occur at 
the respective T=0 thresholds for gas uptake. 

The resulting values of the ground state density, lattice 
constant and chemical potential, pc, dc and p,c , are 
reported in Table I. 

TABLE I. Density, NT lattice constant, and chemical po- 
tential at threshold for gas uptake. 





^He 


Ne 


H2 


D2 


Ar 


CH4 




0.215 


0.327 


0.277 


0.275 


0.273 


0.270 


dc (A) 


17.024 


17.069 


17.166 


17.154 


17.330 


17.373 


(K) 


-380.7 


-927.5 


-480.7 


-621.4 


-1270 


-1290 



Note that the chemical potential at this threshold is 
about 200 K lower than the value —281.5 K found in 
the absence of dilation. This is a measure of the greater 
incentive to absorb within this environment than within 
an undilated NT bundle. Such a dramatic result should 
be testable by thermodynamic measurements, either ad- 
sorption isotherms or direct measurements of the heat of 
adsorption. We note that the computed value of p,c is 
higher than the corresponding threshold, near ^ — 600K, 
for absorption within the tubes (which is relevant only if 
the tubes are open) or in the grooves formed on the exter- 
nal surface of the bundle. The difference is a consequence 
of the much lower zero-point energy in these sites than 
in the IC's. Figure |^ shows the various contributions to 
the energy per particle as a function of p. 

Figure ^ presents our predictions for the dependence of 
d on p for various gases. The isotope dependence arises 
because the smaller mass of H2 implies a larger ZPE. 
Hence, in the D2 case the lattice is dilated less and the 
threshold values of p and the chemical potential are both 
lower than in the H2 case [0. Note that the swelling 
is predicted to be ~ 1% change in the value of d at the 
threshold for gas uptake. This is measurable by X ray 
diffraction from the NT lattice |l^. Table I reports re- 
sults of the effects of dilation on various gases. Two of 
these, He and Ne, fit nearly perfectly in the undilated 
lattice; hence the dilation d is less than 0.5% and the 
increase in binding is small (~ 3% for Ne). For H2 , Ar 
and CH4 , the energetic consequence of the dilation is 
significant. One factor to bear in mind is the sensitivity 
of these results to the potential parameters. As exam- 
ples, in the case of both H2 and CH4, a 2.5 percent in 
the gas-carbon length parameter ( crac) results in about 
a 25 percent increase in the magnitude of pc- 

The nanotubes should exhibit other consequences of 
the H2 confinement. One, possibly surprising, finding is 
a significant shift in the breathing mode frequency of the 
tubes. To evaluate this, we consider the degree of free- 
dom corresponding to the variation in radius of the tubes, 
R. "Breathing" corresponds to a uniform radial expan- 
sion and contraction of a tube. We evaluate the shift by 
adding an extra term 5e{R) = ~ RoY to the previ- 

ous expression for the total energy per unit length, Eq. ^. 
Here i?o is the equilibrium radius in the absence of the H2 
and 7 = 2.5-10^ K A^'^ is the force constant derived from 
the breathing mode frequency of the unperturbed lattice 
|]l9| . Such a change in R yields a change in potential en- 
ergy experienced by the H2 molecules, which we compute 
from the formulation of Stan et al. jl^. As expected by 
analogy with the dilation problem, there ensues a small 
decrease in the tube radius (of order 0.01 A). Much more 
dramatic is the effect on the breathing frequency, which 
can be evaluated by computing the second derivative of 
the total energy with respect to R. Our calculations yield 
a 3% change in the breathing frequency at the threshold 
(p = Pc) for H2 uptake at T = 0. This should be easily 
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measured spectroscopically |T^. The magnitude of the 
shift is a manifestation of the sensitivity of the H2 po- 
tential to the positions of the nearby carbon atoms; it 
is comparable to the shift of this frequency due to tube- 
tube interactions. 

We may also make some predictions concerning the fi- 
nite temperature phase diagram of this unusual system. 
Below a critical temperature Tc there occurs a region of 
two-phase coexistence, between a low density quasi-lD 
"vapor" phase, in a nearly unswoUen NT environment, 
and a high density phase in swollen ICs. At very low 
T, this latter phase is essentially the ground state fluid, 
with a negligible thermal excitation so its chemical po- 
tential satisfies /x(pc, T) ~ niPc, 0). The vapor density at 
coexistence may be determined by equating the chemi- 
cal potentials, a procedure which is particularly simple 
at very low T. Since this vapor phase is dilute, its chem- 
ical potential is well approximated by the ID ideal gas 
relation fi = fc^TlnpA -|- ei(d) where A is the de Broghe 
wavelength of the molecules. By equating these we de- 
termine the curve: p^.X = e'^^^'^''^^°)~'i(''''^' describing 
coexistence of the low density phase and the high density 
phase. At temperatures above a critical value T^, ji^p) 
will no longer exhibit a minimum and only one phase will 
be possible. 

Lacking a better theory, we estimate Tc by mapping 
our system to an Ising system with weak long range forces 
with hamiltonian H = j crif^j, for which the 

mean field approximation is exact and provides a critical 
temperature: T^. = J/Zc^-The mapping is easily found by 
identifying J/N with the effective NT-mediated H2— H2 
interaction energy per particle: Veff = 2[e — ei{do)]. For 
^He, Ne and H2 we find respectively = 5.4 K, 63.8 K 
and 398.4 K. 

To summarize, we have proposed that lattice dilation 
plays an important role in gas uptake. Because the stim- 
ulus to this dilation is the resulting increase in absorbate 
binding energy, the expansion is very sensitive to the gas- 
tube interactions. The inter-tube and intermolecular in- 
teractions also play an important role in determining the 
uptake threshold condition. We have used plausible, but 
uncertain, semiempirical models for these interactions. 
Experimental investigation (thermodynamic, diffraction 
and Raman scattering) will shed light on these assump- 
tions. The existence and nature of this dilation phe- 
nomenon and accompanying transition are robust con- 
clusions, evidently not dependent on the approximations 
used in the interactions. 
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